Dielectric properties of ultrathin metal films around the percolation threshold 
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We report on optical reflection measurements of thin Au films at and around the percolation 
threshold (film thickness 3 to 10 nm) in an extremely broad spectral range from 500 to 35 000 cm -1 
(0.3 - 20 /im). Combining spectroscopic ellipsometry and Fourier-transform infrared spectroscopy, 
the dielectric properties of the films can be described over the whole frequency range by Kramers- 
Kronig consistent effective dielectric functions. The optical conductivity of the films is dominated by 
two contributions: by a Drude-component starting at the percolation threshold in the low frequency 
range and a plasmon in the near-infrared region, which shifts down in frequency with increasing film 
thickness. The interplay of both components leads to a dielectric anomaly in the infrared region 
with a maximum of the dielectric constant at the insulator-to-metal transition. 

PACS numbers: 73.50.-h, 78.20.Ci, 71.30.+h, 73.25. +i 



I. INTRODUCTION 



II. EXPERIMENTAL DETAILS 
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Although for many applications closed metal films as 
thin as possible are desired, the dielectric properties 
of percolating metal films around the insulator-to-metal 
transition are not well understood. Thick continuous 
metal films show a behavior similar to bulk material 1 
and can therefore be well described by the Drude-model 
when corrections for size effects are considered. 2 With 
decreasing thickness, the films become granular, and be- 
low the percolation threshold the metallic behavior dis- 
appears as far as the electrical transport is concerned. 
In principle one can try to simulate semi-continues metal 
films with effective medium (EMA) theories^. How- 
ever it was shown that EMA models fail to describe the 
dielectric properties of ultra-thin metal films 5 '-. In gen- 
eral this transition region was not yet studied in detail 
over a broad frequency range; nevertheless, these inves- 
tigations should eventually provide the data to link the 
abrupt change in dc-conductivity with the observed shift 
of plasmon resonances in the visible. 



On the low frequency side of the electromagnetic spec- 
trum, percolation theories deal with an "idealized" dc- 
conductivity or the electrical behavior at audio and ra- 
dio frequencies up to some MHz£ On the other side of 
the spectrum - in the visible and ultraviolet - a large 
number of investigations have been carried out to un- 
derstand the optical properties of arrays of clusters at 
surfaces. 9 In between, however - in the infrared and far- 
infrared spectral region - little has been done so far. Here 
we report on Fourier-transform infrared reflection spec- 
troscopy and spectroscopic ellipsometry measurements of 
thin Au films on Si/Si02 covering the infrared to ultravi- 
olet spectral range between 500 to 35 000 cm" 1 (60 meV 
to 4.3 eV), corresponding to a wavelength of 280 nm to 
20 /im. The film thickness d was varied between 3 and 
10 nm, i. e. from well below to well above the insulator- 
to-metal transition. 



Thin Au films were prepared by an electron-beam 
heated effusion cell on clean 0.55 mm thick Si(100)- 
substrates covered by a 200 nm Si02 layer. The sub- 
strates were polished on both sides and held at room 
temperature during evaporation. The preparation was 
performed in ultra high vacuum (UHV) at a base pressure 
of 1 • 10~ 7 Pa. The deposition rate was about 1 A/min 
and measured prior to each film-deposition by a quartz 
crystal microbalance. Keeping the gold-flux constant the 
film thickness could be controlled by the evaporation- 
time. The film thicknesses vary from d — 3.4 nm to 
9.5 nm. 

All optical experiments were performed at room 
temperature. The reflectivity measurements in the 
infrared (IR) were carried out by a Bruker IFS 66/s 
Fourier transform infrared (FTIR) spectrometer in the 
range 500-12 000 cm -1 employing a nitrogen cooled 
mercury-cadmium-telluride (MCT) detector. All spectra 
were recorded with a resolution of 0.5 cm -1 and 64 av- 
eraged scans. As reference a thick gold mirror was used. 
To compare measured spectra with model calculations, 
etalon effects (multireflection within the substrate) 
were smoothed out. For the ellipsometry measurements 
a Woollam variable angle spectroscopic ellipsometer 
(VASE) was utilized. The experiments were carried out 
in the spectral range between 6000 and 35 000 cm -1 
(0.75 to 4.3 eV, corresponding to a wavelength of 280 nm 
to 1.7 /im) with a resolution of 170 cm^ 1 and an angle 
of incident varied in steps of 5° from 35° to 75°. 



III. RESULTS AND ANALYSIS 

The thickness dependent morphology of thin Au-films 
on Si/Si02 is well investigated^. In Fig. Q] atomic force 
microscopy (AFM) images at the same scale from a for- 
mer investigation of 3 nm, 5 nm, and 7 nm thick gold 
films are presented 6 . 
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FIG. 1: (Color online) AFM images of 3 nm, 5 nm, and 7 nm 
thick gold films on Si/Si02- All three images are presented 
at the same scale. With increasing film thickness the cluster 
size increases and the films become smoother. 

With increasing film-thickness the size of the islands 
increases, they get closer to each other and the films 
become smoother. The root-mean-square roughness of 
the 3 nm film is about 1 nm, and the average clus- 
ter size is about 4 nm. The percolation threshold it- 
self can not be determined from the morphology. Al- 
though it is important to know the morphology of the 
films, one has to stress the fact that percolation theories 
deal with scaling laws and therefore the overall electri- 
cal and optical behavior should be independent of the 
specific morphology^. 

Prior to the film preparation, the bare Si/Si02 sub- 
strates were characterized by ellipsometry and the ob- 
tained optical parameters then kept constant in the fur- 
ther modeling of the films. As an example, in Fig. [5] 
the measured ellipsometric angle \I/ for a 8.1 nm thick 
Au film is shown as a function of frequency and an- 
gle of incidence. Several attempts were made to model 
the obtained ellipsometric angles: It was not possible, 
for instance, to model the granular films by an effective 
medium approximation (EM A), even in the case where 
the volume fraction was gradually decreased over the film 
thickness. Taking the granular Au layer as a homoge- 
neous effective layer with fixed thickness, it was possible 
to get reasonable agreement to the experimental data 
by a point-by-point fit with MSE-values around 6; but 
it turned out that the received dielectric functions were 
not Kramers-Kronig consistent. 

To solve this problem, we combined the results ob- 
tained by ellipsometry and IR-reflectivity. Both sets 
of data were simultaneously analyzed with the p rog ram 
package RefFIT. For details on this tool see Ref. [lCj. By 
means of a variational dielectric fit within the program, 
Kramers-Kronig consistent dielectric functions could be 
obtained, which perfectly reproduces both the ellipso- 
metric angles (continuous line in Fig. [2]) as well as the 
IR-reflectivity. In Fig. [31 for example, the measured re- 
flectivity in the IR-region (dashed line) of a 6.5 nm thick 
Au film with its corresponding fit (solid line) is shown to- 
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FIG. 2: (Color online) Measured ellipsometric angle ^ for 
a 8.1 nm thick Au film on a Si/Si02 substrate for several 
angles of incidence (dashed lines). The solid lines are the 
corresponding model fits. For details see text. 
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FIG. 3: (Color online) Measured IR-reflectivity of a 6.5 nm 
Au film (dashed line) together with the reflectivity calculated 
from the Kramers-Kronig consistent dielectric functions ob- 
tained by a simultaneous analysis of the IR- and ellipsometric 
data. The inset shows the assumed layer structure. 

gether with an extrapolation to the ultraviolet obtained 
from the ellipsometric measurement. The dominating os- 
cillation in the reflectivity is due to interferences caused 
by the 200 nm thick SiC>2 layer. 

The real part ei(w) of the dielectric function as received 
from the above analysis is shown in Fig.Uas a function of 
frequency for varying film thickness. From the imaginary 
part 62(10) the optical conductivity^ was calculated by: 

01 (u>) = ue e 2 (uj), (1) 

with eo the permittivity of free space. It is displayed in 
Fig. [5] in a double logarithmic fashion. 

The dielectric functions <j\(uj) and ei(w) are described 
by a superposition of Drude- and Lorentz-terms*ii The 
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FIG. 4: (Color online) The effective dielectric constant of 
the gold films. Note the non-monotonous behavior at low 
frequencies as the film thickness d increases (indicated by the 
arrow). From a gold-film thickness of 8.f nm on ei(w — » 0) is 
negative, indicating the metallic behavior of the films. 
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FIG. 5: (Color online) Optical conductivity of gold films of 
different thickness as indicated. The percolation occurs be- 
tween 6.5 and 7.0 nm where the sign in the slope of ai(u) 
changes at low frequencies. 



parameters obtained from the fits for the different film 
thicknesses are listed in Table |U 



IV. DISCUSSION 

A. Percolation threshold 

From the effective conductivity <7i(w) as well as from 
the effective dielectric function ei(ui) of the gold layers, 
the different films can be clearly divided in two regimes: 
The continuous films in the thickness range d > 7 nm 
and the granular films with d < 6.0 nm. The percolation 
threshold d c falls right between these two values. 



Well above the percolation threshold, the conductiv- 
ity decreases with frequency and the static permittivity 
ei(u> —> 0) is negative, both indicating a typical metal- 
lic behavior. The low-frequency behavior of the metal- 
lic films can be fitted by a simple Drude modelii with 
plasma frequencies lo p and scattering rates 7 = 1/ (2ttct) 
as summarized in Table |H In addition, the fit parameters 
of two Lorentz oscillators (discussed below) are listed. 
Naturally, the reflectivity of the evaporated films is lower 
compared to bulk Au, because the effective electron den- 
sity is reduced and the scattering rate 7 becomes larger 
with decreasing film thickness. While the former effect 
comes from a surface dipole laye r 12 i 13 , the latter is a well- 
known consequences of the classical size effect. Both ef- 
fects become more pronounced as d is reduced since the 
surface-to- volume fraction is enhanced. 14 ! 15 ' 16 

Below the metal-to-insulator transition (MIT) the 
static conductivity a\{uj — > 0) vanishes and ei(w) is pos- 
itive. The optical behavior in this regime can be ex- 
plained by the absence of the Drude contribution and 
two resonance peaks in the mid-infrared. Below the MIT 
we picture metallic islands which interact capacitively. 
While the dc-conductivity is zero, the optical conductiv- 
ity depends on the capacitive coupling of the islands and 
thus increases with frequency From the AFM images of 
comparable samples (see Fig. [T]) and other very exten- 
sive studies^ of the morphology of Au on Si/Si02, the 
thickness dependence can be explained as follows: the 
metallic clusters are separated from each other forming a 
condenser. With increasing effective film thickness d the 
islands are closer to each other and also larger. While 
their surface increases, the spacing decreases and tends 
to zero at the percolation threshold. This leads to a rise 
of the capacitive coupling and an enhanced optical con- 
ductivity (Fig. [5]). For a detailed discussion from a more 
theoretical point of view see e. g. Ref. Ii"7lll8l . 

B. Dielectric constant 

As far as the thickness dependence of the low-frequency 
permittivity is concerned, one would expect that e±, 
which is positive for d = 0, gradually decreases with d and 
becomes negative at the MIT. However, a distinctively 
different behavior is observed in our experiments. As in- 
dicated by the arrow in Fig.[2J ei(w) first rises with film 
thickness d, goes through a maximum at some critical 
concentration d c which is assumed to be around 6.0 nm, 
and then decreases rapidly; e± becomes negative only at 
considerably larger d, above 7 nm. 

In principle, this behavior is known from the low- 
frequency conductivity of percolating networks. Efros 
and Shklovskii investigated such systems theoretically 
and predicted a divergence of the static dielectric con- 
stant 

e 1 (0,d)oc(d c -d)- s (2) 
at a critical thickness d c J£ For any small but finite fre- 
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TABLE I: The Drude- and Lorentz-parameters obtained from the fit to the model. The plasma frequency is denoted by uj p and 
7 = l/(27r ct) is the scattering rate with c the speed of light and r the relaxation time. The bulk values are taken from Ref. QJ. 
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to a two-dimensional system since the data can also be 
explained with the same critical thickness but the three- 
dimensional critical exponent s — 1 (not shown). The 
dashed lines indicate the maximum value of the dielectric 
constant ei(1000 cm -1 , d c ) for two and three dimensions 
according to Eq. 

For a more accurate analysis directly at the percolation 
threshold and a quantitative prove of equation ad- 
ditional data with closer steps in thickness are required. 
It is interesting to note that the 7 nm film still exhibits 
a positive e\ (cf. Fig. 2]) albeit it possesses a Drude com- 
ponent. 



Plasmons 



FIG. 6: (Color online) Divergence of the dielectric constant at 
1000 cm- 1 . The solid line is a fit to the data: e\ oc (6 — d)~ 1,3 
for d < 6 nm and linear above. The dashed lines represent 
the predicted maximum after Eq. ([3} for two and three di- 
mensions. 

quency to, the divergence simply becomes a maximum: 

( \ lu 
— — I > ( 3 ) 
u>e e s J 

where a m is the real part of the conductivity of the metal- 
lic fraction and e s the real part of dielectric constant of 
the substrate. For a two-dimensional system the critical 
exponents in Eqs. ([2]) and ([3]) are s — 1.3 and u = 0.5. 
For three dimensions s = 1 and u = 0.62. 18 Figure [S] ex- 
hibits the dielectric constant as a function of film thick- 
ness measured at a frequency of 1000 cm -1 ; it represents 
more or less the static value since ei is almost frequency 
independent below this value. The measurements fit to 
the expected behavior with a maximum at the critical 
thickness d c somewhere around 6 nm. The zero-crossing 
of ei(u) — 0, d) occurs at do slightly above 7 nm, where 
the Drude-component in the films start to develop (cf. 
Table fl]). The solid line corresponds to Eq. with 
d c = 6 nm and s = 1.3 but does not unambiguously refer 



For the description of the behavior at higher frequen- 
cies additional Lorentz-oscillators have to be considered. 
To illustrate how the conductivity is built up by the dif- 
ferent components, the conductivity of the 9 nm film is 
shown in Fig. [7] together with its different contributions. 
With increasing film thickness the two oscillators shift to 
lower frequencies, become weaker and narrower (cf. Ta- 
ble U). 

In the following we will first concentrate on the sec- 
ond oscillator at higher frequencies which is caused by 
the ensemble of single particles. It is a property of the 
ensemble and has to be distinguished from free electron 
oscillations in single metal clusters. The incident elec- 
tric field is modified by the polarizability of these parti- 
cles which - most important - interact ] 19 ' 20 In literature 
this transverse-mode plasmon is referred to as "Maxwell- 
Garnett resonance" or "optical conduction resonance". 
It shifts to lower frequencies upon stronger interaction as 
the metal-islands become closer and denser packed. As- 
suming spherical clusters with a diameter much smaller 
than the used infrared wavelength, the position of this 
plasma absorption can be described by the following con- 
ditioni 9 -: 

^(u) = ei = ■ (4) 
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FIG. 7: (Color online) Decomposition of the measured con- 
ductivity of the 9 nm film in its Drude- and Lorentz- 
components as listed in Table [I] 



In this condition, Q is the area fraction of the substrate 
covered by metal, n s is the refractive index of the sub- 
strate and e\ ulk (cu) is the frequency dependent dielectric 
constant of the metallic fraction. For a given Q, we can 
calculate the right-hand side of the equation and obtain 
a certain dielectric constant e\. The frequency at which 
the obtained resonance condition e\ ulk (ui) = t\ is fulfilled 
gives the position at which the plasmon appears. 

If we use the frequency dependent values e\ ulk (oj) of 
bulk gold as calculated from^L and set n s — 1.45 for 
the SiC>2-substrate, we get the solid line in Fig. [5J We 
can now compare this line with the measured plasmon 
frequencies in dependence of film thickness as listed in 
Table |T] as second Lorentz-oscillator ujq . Here we assume 
a linear relation between the film thickness d and the area 
fraction Q covered by the metal. As seen from Fig. [8j 
the plasmon shifts to lower frequencies with increasing 
film thickness d, nicely following the behavior given by 
Eq. Q. 

A closer inspection yields small deviations of the mea- 
sured plasmon frequencies to lower values, which we as- 
cribe to the influence of the substrate. The optical prop- 
erties of the 200 nm oxide layer might be influenced by 
the underlaying silicon, which has a much higher refrac- 
tive index. Taking a somewhat larger n s in Eq. ^ shifts 
the calculated values to lower frequencies, i. e. closer to 
the experimental data. 

As can be seen from Fig. [5] for Q — > 1 the calculated 
plasmon frequency shifts to zero and thereby mutates to 
a Drude-peak. It is known that the Maxwell-Garnett the- 
ory breaks down above the percolation threshold, which 
in our case is at about Q = 0.6, therefore the extrap- 
olation to Q — > 1 is disputable. Nevertheless, the ex- 
perimentally found position of this plasmon can still be 
described quite well within this theory even above the 
percolation threshold (see Fig [5}. 

Recently de Vries et al£2- interpreted their ellipsomet- 
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FIG. 8: (Color online) The position of the plasmon shifts to 
lower frequencies with increasing film thickness. The solid 
line is calculated from Eq. (j4| with the values of bulk gold 21 
and belongs to the Q-axis. The dots show the experimental 
position of the plasmon for films with thicknesses as given 
by the upper-axis. The proportional relation between the 
two abscissas is arranged in the way that the slope of the 
experimental and calculated values fits best. 



ric measurements on thin silver films in the visible and 
near-infrared spectrum exactly in this way: the resonance 
frequency of a localized plasmon shifts to lower frequen- 
cies with increasing film thickness becoming zero at the 
percolation threshold. There the relaxation time exhibits 
an abrupt increase indicating the transition to a macro- 
scopic conducting state. 

Our measurements, which extent to the far infrared, 
reveal a different picture. At the percolation threshold 
we can clearly see both the development of a Drude-peak 
and the presence of the plasmon. Even well above the 
percolation threshold, the position of this plasmon can 
still be described satisfactorily by Eq. Q, but it gradu- 
ally dies out with increasing film thickness d. This can 
in principle be interpreted in the way that the rough sur- 
face of the percolated films still show dipole interaction 
or that dielectric inclusions in the film start to interact 
with each other, as it has been considered by Cohen et 
al.—. In a direct comparison of the conductivity of the 
calculated and measured plasmon the Lorentzian-line of 
the latter one is broadened (data are not shown). This 
comes from the fact that the particle size is smaller than 
the mean free path of the conduction electrons. As seen 
before, this size effect also leads to a large scattering 
rate in the Drude component of the percolated films (cf. 
Table [TJ) . Therefore the Maxwell-Garnett theory gives 
correctly the resonance position but the calculated plas- 
mon band is wider than expected. The occurrence of this 
effect for thin films was already predicted by Doremus 19 . 

As discussed above the plasmon referred to as second 
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FIG. 9: Spectral weight and electron density N e of the Drude- 
and the two Lorentz-oscillators as function of the film thick- 
ness. The dashed line represents the electron density as de- 
termined for bulk gold. 



frequencies is now a direct consequence of the two com- 
peting contributions: below the percolation threshold the 
shift of the plasmons to lower frequencies with increasing 
film thickness leads to an increasing ei(cj). At the per- 
colation threshold d c the Drude-peak starts to develop, 
which adds a strong negative component to the dielectric 
response leading to a maximum in the dielectric func- 
tion. Eventually t\ changes sign at do which is slightly 
above d c . This interpretation is supported by the anal- 
ysis of the spectral weight and the corresponding elec- 
tron density N e shown in Fig. [5] The spectral weight 
increases linearly with film thickness reaching the bulk 
value at about 10 nm. No indication of the MIT be- 
tween 6.0 nm and 7 nm can be identified, i. e. the abrupt 
change in dc-conductivity is not reflected at higher fre- 
quencies. There a monotonous transfer of spectral weight 
from the plasmons to the Drude-peak is observed. The 
linear increase of the electron density simply describes 
the growing amount of metal. A similar result for thin 
silver films was observed by de Vries et al^. 

V. CONCLUSION 



Lorentz oscillator in Table U is a property of an ensem- 
ble of single clusters. At low coverage only this peak 
is seen. At higher coverage a second peak at lower fre- 
quencies appears. This mode is an indication for the 
formation of aggregates^. With increasing film thick- 
ness both plasmon peaks shift to lower frequencies, they 
broaden and the splitting between them increases. In 
the classification given by Kreibig et al£z the sample un- 
dergoes a transition form category II (separated single 
clusters in full statistical disorder) to category IV (vari- 
ous kinds of aggregates, various next neighbor distances, 
plus single particles). At the percolation threshold, when 
additionally the Drude component appears, they can not 
be clearly separated anymore (see Fig. [7|). 

The maximum in ei(w) at the critical thickness at low 



Combining FTIR spectroscopy and spectroscopic ellip- 
sometry the effective dielectric function of thin Au films 
around the percolation threshold could be obtained over 
a very broad frequency range from the infrared up to 
the UV. The optical properties of the films can in prin- 
ciple be described by two contributions: a plasmon in 
the near-infrared, which shifts down with increasing film 
thickness and than slowly dies out above the MIT, and 
a Drude peak, which starts to develop at the MIT and 
than rapidly increases with film thickness. The interplay 
of both components leads to a dielectric anomaly, known 
from percolation theory: from dc up to a few thousand 
wavenumbers £1(0-") exhibits a pronounced maximum at 
some critical thickness. 
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